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Abstract. We present an asymptotic expansion for quaternionic self- 
adjoint matrix integrals. The Feynman diagrams appearing in the expan- 
sion are ordinary ribbon graphs and their non-orientable counterparts. 
The result exhibits a striking duality between quaternionic self-adjoint 
and real symmetric matrix integrals. The asymptotic expansions of these 
integrals are given in terms of summations over topologies of compact 
surfaces, both orientable and non-orientable, for all genera and an ar- 
bitrary positive number of marked points on them. We show that the 
Gaussian Orthogonal Ensemble (GOE) and Gaussian Symplectic En- 
semble (GSE) have exactly the same graphical expansion term by term 
(when appropriately normalized), except that the contributions from non- 
orientable surfaces with odd Euler characteristic carry the opposite sign. 
As an application, we give a new topological proof of the known dual- 
ity for correlations of characteristic polynomials. Indeed, we show that 
this duality is equivalent to Poincare duality of graphs drawn on a com- 
pact surface. Another application of our graphical expansion formula is 
a simple and simultaneous (re)derivation of the Central Limit Theorem 
for GOE, GUE (Gaussian Unitary Ensemble) and GSE: The three cases 
have exactly the same graphical limiting formula except for an overall 
constant that represents the type of the ensemble. 
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1. Introduction 

The purpose of this paper is to establish an asymptotic expansion for 
quaternionic matrix integrals in terms of Feynman diagrams which exhibits 
a duality between quaternionic self-adjoint and real symmetric matrix inte- 
grals. 

Recent developments in the theory of random matrices exhibit particu- 
larly rich structures. Although originally introduced by Wigner as a model 
for heavy nuclei, random matrices appear almost ubiquitously in mathemat- 
ics and physics. Applications pertain, for example, to number theory [11], 
combinatorics [1], probability theory [20], algebraic geometry [12] and 
quantum gravity [22]. 't Hooft's discovery that quantum chromodynam- 
ics (QCD) simplifies in the limit where the number of colors N (i.e. gauge 
group SU(N)) is large relied on a graphical expansion in terms of "fat" 
or "ribbon" graphs [19]. While the three major ensembles (Gaussian Or- 
thogonal, Gaussian Unitary, and Gaussian Symplectic) have important ap- 
plications, simple graphical expansions in terms of ribbon graphs and their 
non-orientable counterparts are currently available only for the GUE and 
GOE models. What type of graphs should be used for the GSE models? 
The non-commutativity of matrix entries has defied earlier attempts to find 
an efficient asymptotic expansion formula. This article provides the missing 
graphical expansion for symplectic matrix integrals. 
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A 2N x 2N symplectic matrix M 

(1.1) MCM T = C 

(C = -C T = C~ f ), may be written as 1 

(1.2) M = I 2x2 ® S + %a x <g> A 1 + icr 2 <g> A 2 + m 3 ® A 3 , 

where S* is a symmetric and A±, A 2 , A 3 are antisymmetric N x N matrices. 
In turn, representing (ia x ,ia 2 , ia 3 ) as quaternionic units (i, j, k) so that 

(1.3) M = S + iA 1 +jA 2 + kA 3 , 

the problem is to write the Feynman rules for a quaternionic matrix integral. 

't Hooft's insightful fat graph notation arises from identifying each index 
a and b of an iV x N matrix (X ab ) as an edge of a ribbon. The most naive 
approach therefore, is to append additional indices a,/3 = 1,2 correspond- 
ing to the 2x2 Pauli matrices in (1 .3). Unfortunately this method leads only 
to partial results [9]. The next most obvious idea is to employ the quater- 
nionic realization, in much the same way that an integration over hermitian 
matrices X = X^ is immensely simplified when one views the entries as 
complex numbers rather than X = S + iA (where the real matrices S and 
A are symmetric and antisymmetric, respectively). The drawback is that 
the quaternions do not commute, so a non-commutative Feynman calculus 
is required. Indeed it is possible to generalize the usual Schwinger trick 
to quaternionic source terms, and represent a Gaussian symplectic integral 
as non-commutative quaternionic differentiations. The result is a sum over 
both orientable and non-orientable ribbon graphs, and is easily verified to 
agree with ours for simple graphs. This method is described in Appendix A. 

Our result for the quaternionic matrix integral is 



(1.4) log 



/ J[dX] exp (- N tr X 2 + ^ ^ tr 
\ J[dX}exp(-NtrX 2 ^) 



11* 



ki i Aut < r >i i ' 

Exact conventions are given later, at present it suffices to indicate that the 
sum is over all graphs T drawn on compact orientable and non-orientable 
surfaces Sr, and x(Sr) is the Euler characteristic of the surface SV uniquely 
defined by the graph T. Our proof of this result is based on viewing the 

1 The Pauli matrices are 

0l\ / -i \ / 1 



ai = l 1 ) ' a2 = { i ' CT3 = I 1 
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GOE integral over real symmetric matrices as fundamental. The crucial 
observation is that the contribution of any given graph T in the GOE, GUE 
and GSE is a topological invariant of the surface Sr with fr marked points 
on it, where fr denotes the number of faces of the cell-decomposition of 
S r defined by the graph T. The connectivity of the space of triangulations 
of two dimensional surfaces then allows any graphical contribution to be 
calculated from a simple representative graph for any given topology. 

Writing the results for all three ensembles in a uniform notation (see (5.2)) 
makes a new duality 

GOE < — ► GSE 



(1.5) GUE < — ► GUE 

GSE < — ► GOE , 

manifest. The middle line for the GUE is a (trivial) self-duality. The tilde 
on the right hand side indicates that equality holds upon doubling/halving 
the matrix size and an overall sign change for contributions of graphs where 
the Euler characteristic x(^r) is °dd. 

Asides from the computational utility of this duality, as an immediate 
corollary, it provides a simple proof of the known duality for A; -fold corre- 
lations of characteristic polynomials of N x N matrices [2, 4, 5, 13]. For 
the GUE model expressed in terms of ribbon graphs, this N-k duality [4, 5] 
is precisely Poincare duality of graphs drawn on a compact oriented sur- 
face. Similarly, the relation between GOE and GSE correlations stems from 
the combination of Poincare duality, this time including non-orientable sur- 
faces, and our new graphical expansion formula. It is interesting to note that 
the machinery of fermionic integrations employed in [4, 5] is equivalent to 
a very simple switch from a graph on a surface to its dual graph. 

If we reduce our integral (1.4) to a symplectic Penner model by setting 
t\ = h = and 

tj = -zi~ 1 , j>3, 
then we can explicitly compute the asymptotic expansion in z of 



(1.6) lim log 



l RN \m\ 2a uLe W [~ k 4)dk t 



utilizing the Selberg integration formula and the asymptotic analysis tech- 
nique of [14], where a is either a positive integer or its reciprocal. We 
demonstrate that the duality we found for GOE and GSE of (1.4) extends 
to the same type of duality between an arbitrary positive integer a and l/a 
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for (1.6) with the sign change for all terms with odd powers of iV in the 
asymptotic expansion. 

The orthogonal Penner model gives the orbifold Euler characteristic of 
the moduli spaces of smooth real algebraic curves with an arbitrary num- 
ber of marked points [6]. Also, the original Penner model [16] provides 
the orbifold Euler characteristic of the moduli spaces of pointed algebraic 
curves over C [7]. Therefore it is natural to ask what the symplectic Penner 
model gives. Our duality shows that actually the symplectic Penner model 
is identical to the orthogonal Penner model, except for a change of the ma- 
trix size and an overall sign for contributions from surfaces of odd Euler 
characteristic. This equivalence has been observed independently in [15]. 

Our final application, is a simple and unified (re)derivation of the Cen- 
tral Limit Theorem for Gaussian random matrix ensembles. This result fol- 
lows as a direct consequence of 't Hooft's original large N limit in which 
planar ribbon graphs dominate: we derive a precise limiting formula for 
GOE, GUE and GSE matrix ensembles in terms of planar two-vertex rib- 
bon graphs. The formula is the same for all three ensembles, except for an 
overall constant. 

The material is organized as follows: In Section 2 we introduce the ma- 
trix integrals studied in this paper. Our conventions for the topological data 
of surfaces are given in Section 3 as well as our theorem and its proof for 
the graphical expansion of matrix integrals. Examples, including a compar- 
ison with the first few terms of the Penner model, are given in Section 4. 
The new duality appears in Section 5. The central formula of this paper is 
equation (5.3) which gives the graphical expansion for the GOE, GUE and 
GSE simultaneously in a manifestly duality invariant form. Its application 
to characteristic polynomial duality is in Section 6. The extended version 
of our duality for Penner type models is found in Section 7 while detailed 
derivations of the formulae there are presented in Appendix B. Section 8 
concerns the the Central Limit Theorem for Gaussian random matrix en- 
sembles. In the Conclusions (Section 9) we discuss possible further gener- 
alizations. In particular, the construction of a graphical topological invariant 
of surfaces with marked points necessary for the proof of our main result, 
is rather general and may be applied to higher algebraic structures. 
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2. Matrix Integrals 
The object of our study is the integral over self-adjoint matrices 2 

J[rfX] (/3) exp(-ItrX 2 + Er=i|trX^) 



(2.1) Z ( P\t,N) 



J[rfX] (/3) exp(-itrX 2 ) 



as a function of the "coupling constants" t — (ti, t 2: t 3: . . . ) and the size N 
of the matrix variable X. Here 

(2.2) X = S + Y,*Ai 

i=l 

is built from real, N x N, symmetric and antisymmetric matrices S and 
Ai, respectively. The parameter f3 takes values 1, 2 or 4 depending whether 
we study real, complex or quaternionic self-adjoint matrices and in turn 
Gaussian orthogonal, unitary or symplectic ensembles (GOE, GUE, GSE). 
The imaginary units e% are then drawn from one of three sets 



(2.3) e t e { 



0, P=l, 

{i:i 2 = -l}, [3 = 2, 

_ k:i 2 = f = k 2 = ijk = -1} , /3 = A. 

The self-adjoint condition 

(2.4) Xt = X T = X , el = -a , 

is implied by antisymmetry of the matrices A^. Finally, the measure [dX]^ 
is the translation invariant Lebesgue measure of the vector space of real 
dimension ^N([3(N — 1) + 2) spanned by independent matrix elements of S 
and A,^ This measure is invariant, respectively, under orthogonal, unitary 
and symplectic transformations 

(2.5) X i — ► U^XU, 

where iPU = 1 and U = U + Y,i=i e iUi for real N x N matrices U 
and Ui. 

The matrix integral (2.1) is a holomorphic function in (t± : t 2 , • • • , hm) if 
we fix m > and restrict the coupling constants to satisfy 

Re(t 2 m) < and t 2m +l = ^2m+2 = hrn+3 = • • • = 0. 



2 We employ various normalizations throughout the paper, so it is convenient to divide 
through by the free matrix integral. 



DUALITY FOR MATRIX INTEGRALS 



7 



Under this restriction, Z^(t, N) has a unique Taylor expansion in (ti, t 2 , 
■ ■ ■ , tim-i) and an asymptotic expansion in t 2m as t 2m — > while keeping 
Re(t2m) < 0. Let us introduce a weighted degree of the coupling constants 
by deg(t n ) = n. Then the asymptotic expansion of the truncated integral 
Z^(t, N) has a well-defined limit as m — > oo in the ring 

of formal power series in infinitely many variables with coefficients in the 
polynomial ring of N with rational coefficients [14]. The subject of our 
study in what follows is this asymptotic expansion of Z^{t, N) as a function 
in infinitely many variables. 

3. Graphical Expansion 

The graphs appearing in our asymptotic expansion of the matrix inte- 
grals (2.1) are those drawn on orientable as well as non-orientable surfaces. 
To avoid confusion with an already well-established convention that rib- 
bon graphs are drawn on orientable surfaces, we propose the terminology 
Mobius graphs. Let us recall that a ribbon graph T is a graph with a cyclic 
order chosen at each vertex for half-edges adjacent to it. Equivalently, it is 
a graph drawn on a compact oriented surface S giving a cell-decomposition 
of it. The complement S \ V of the graph r on S is the disjoint union of 
fr open disks (or faces) of the surface. Since a ribbon graph T defines a 
unique oriented surface on which it is drawn as the 1 -skeleton of a cell- 
decomposition, we denote the surface by Sr. 

Similarly, a Mobius graph is drawn on a compact surface, orientable or 
non-orientable, giving a cell-decomposition of the surface. It can be viewed 
as a ribbon graph with twisted edges. A Mobius graph T also uniquely 
defines the surface Sr in which it is embedded. 

Let (5 be the set of connected Mobius graphs. A graph T E (5 consists of 

(i) 

a finite number of vertices and edges. Let v T . denote the number of j-valent 
vertices of T. Then the number of vertices and edges are given by 

(3-D = and e r = \ X>£°. 

j j 

The unique compact surface 5V has / r faces and its Euler characteristic is 
(3.2) X (Sr) = vr - e r + fr • 

We will also need the number of faces with a given number of edges, so de- 
note the number of j-gons in the cell-decomposition of Sr by f^ whereby 
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the Poincare dual formulae to the equations (3.1) are 
(3.3) /r = X)/r ) and e r = \ 



2 

3 3 

A Mobius graph T also determines the orientability of S r and we define 

{1 SV orientable , 
— 1 S r non-orientable . 

By genus of Sr we mean 

(3.5) g(S T ) = 1 - 2"^ X (Sr) ■ 

Thus x(Sr) = 2 — 2g(Sr) for an orientable surface and x(Sr) = 1 — g(Sr) 
if it is non-orientable. We also define the parity of X (Sr) by 

(3.6) = (-I)**) , 

while many results can be written compactly in terms of 

T orientable , 



(3.7) S r = l - (l + (tr) - 



Our main result is: 



1 T non-orientable, X (Sr) odd , 

2 T non-orientable, X (S T ) even • 



Theorem 3.1. The logarithm of the asymptotic expansion of the matrix in- 
tegral Z^> (t, N) is expressed as a sum over connected Mobius graphs: 

(3.8) log(z^,A0) 

= v (-4 + 6/3- (3 2 ) 1 -^ xiSr) (2 - pFp^N* n ,0) 

h ' Aut ( r )i A / j 

€{Q[N])[[ti,t2,t 3 ,...]]. 

Remark: 

1. We define (2 - /3) Sr = i when (3 = 2 and S r = 0. 

i-r v U) 

2. For every connected Mobius graph T, the monomial [[j tj T is a finite 
product of total degree 2e r . 

3. The only reason to consider logZ /3 (t, N) is because it yields a com- 
pact formula (3.8): Z@(t, N) itself has an expansion in terms of graphs 
although a given summand may have a mixture of orientable and non- 
orientable connected components. 
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The automorphism group Aut(r) of a Mobius graph T is a group of 
automorphisms of the cellular complex Sr consisting of vr vertices, er 
edges and fr faces. When Sr is orientable, the group Aut(T) may con- 
tain orientation-reversing automorphisms as well. We note that a cyclic 
rotation of half-edges around a vertex corresponds to the invariance of the 
trace under a cyclic permutation 

(3.9) tr (.\/ ; .\ /,.!/, • • • M n ) = tr (M n M 1 M 2 ■ ■ ■ M n _x), 

and an orientation-reversing flip of vertex with adjacent half-edges corre- 
sponds to the invariance the trace of symmetric matrices 

(3.10) tr (SiSzS 3 ■■ - S n ) =ti (S n S n -iS n -2 ■ ■ ■ Si) 
under reversing the order of multiplication. 



Figure 3.1. Two equivalent Mobius graphs consisting of 
two vertices, three edges, and one face. The graphs are in- 
terchanged by a vertex flip. 

The proof of the theorem involves two main ingredients. The first is to 
view GUE and GSE matrix integrals as the coupling of a singlet or triplet of 
skew-symmetric matrix integrals to the fundamental GOE integral. When 
P = 1, equation (3.8) is the Mobius graphical expansion of a symmetric 
matrix integral [6, 21, 18, 10]: 

an, log^t^^^^n^. 

This formula follows immediately from the fact that the Wick contraction 
of any pair of symmetric matrices S = (Sab) obeys 



(3.12) 



(SabScd) — 8ac8bd + &ad&bc , 
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which is denoted graphically as an edge of a Mobius graph (see Figure 3.2). 
These edges connect vertices of the type 

t- t- N 

(3-13) -^7 tr S J = -^r S aia2 S a2aa ■ ■ ■ S ajai . 

J J ai,...,aj=l 

as depicted in Figure 3.3. The factor (2j) _1 in (3.13) is precisely the one re- 
quired to cancel the over-counting implied by the identities (3.9) and (3.10). 
Therefore the overall weight of any given graph is as quoted in (3.1 1). 



b 
a 



^ab 



+ 



b 
a 



c 
d 



c 
d 



^cd 



Figure 3.2. Propagator for symmetric matrices. 





Figure 3.3. Two vertices connected at the half edge la- 
beled by aia 2 to the half edge labeled by b 2 bi. For the GOE 
a second graph with a twisted edge connecting aia 2 to bib 2 
is also present since the ribbon edges are not directed. 



In contrast, for a pair of antisymmetric matrices A = (A a b) the Wick 
contraction yields 

(3.14) (A ab A cd ) = 5 ac 5 bd - 5 ad 5 bc . 

This is depicted in Figure 3.4. 
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b 
a 



A 



ab 



b 
a 



c 
d 



A 



cd 



c 
d 



Figure 3.4. Propagator for antisymmetric matrices 

The minus sign for untwisted edges will play a crucial role in what fol- 
lows. Notice that since the exponent of the Gaussian part of (2.1) is 

, N (3-1 

(3.i5) -\J2 s i-\ £ (^+5>l), 

a=l l<a<6<7V 1=1 

there are no non- vanishing Wick contractions between symmetric and anti- 
symmetric matrices. Therefore the model (2.1) is a sum over graphs with 
edges of either of the two types (3.12) or (3. 14). Furthermore (AiAj) is only 
non- vanishing when i = j so edges "emitted" at vertices are only connected 
when they carry the same imaginary units ej. 

If we call e = 1 and A = S, a j-valent vertex now looks like 

(3.16) 

h * * = £ e « A « = £ PK 2 V' Qi) tr II A ak e ak . 

«=0 {(ai,...,a 3 )}/~ ^ fe=l 

The last sum runs over all j-tuples of integers 0, 1, — 1 modulo the 
equivalence relation 

(3.17) 

(cti, . . . , aj) ~ ((3i, . . . , /%) iff tr (A ai • • • A a .) = tr (A^ • • • A 0j ) . 
The multiplicities p(ai, ... , a-,) are non-zero or/j when the product 

(3.18) e ai ■ '-e aj = ±1 . 

In other words, the vertex (3.16) is real. Observe that the numbers 

(3.19) ^ 

p{ai, . . . ,OLj) 

count the number of automorphisms of a vertex with a given configuration 
of units e a sprinkled at every edge, with respect to rotations and vertex flips. 
For example, for j3 — 4 and j = 3 

itrX 3 = itr(l.S) 3 
6 6 K J 
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+ X - tr (l.S- i.A x i.Ai + l.S j.A 2 j.A 2 + 1.5 k.A 3 k.A^j 
(3.20) + i tr (i.^ j.A 2 fc.A 3 ) . 

On the first line, the trace of three identical symmetric matrices can be cy- 
cled and reversed, so 2j = 6 is the correct automorphism factor. On the 
second line the terms with a single symmetric and a pair of antisymmetric 
matrices can only be flipped yielding two automorphisms. The term on the 
third line has no automorphisms. We will keep track of the matrix type by 
"sprinkling" units {1, e^} over the set of Mobius graphs (see Figure 3.5). 

Orchestrating the above observations yields the following: 

Lemma 3.2. The matrix integral log JV)) may be computed as a 

sum over Mobius graphs Y with weight 

(3 - 21) fern 

multiplied by a single power oftj for each j-valent vertex in T. The factor 
fir is calculated by 

1. Writing down all possible configurations of units e a G {1, ei, . . . e /3 _ 1 } 
at each vertex such that their product is ±1. 

2. Counting these signed configurations with an additional minus sign 
for every twisted edge with imaginary units e$ at each end. 

The proof of this lemma follows from our previous remarks and by noting 
that including e a from the vertices, the Wick contractions become (no sum 
on i) 

(3.22) (l.S ab l.S cd ) = S ac S bd + 5 ad 5 bc , 

(3.23) (ei.A iah ei.A icd ) = -5 ac 5 bd + 5 ad 5 bc . 

A sample computation of fir is given in Figure 3.5. 

Our proof is completed by computing the numbers fir, which requires a 
second major ingredient: 

Lemma 3.3. The quantity fir is a topological invariant of the n-punctured 
surface Sr \ {pi,P2, ■ ■ ■ ,Pn}> where the removed points are all distinct and 
n = fr is the number of faces of the cell decomposition defined by T. 

The proof of this crucial Lemma is almost a triviality. Graphs are iden- 
tified if they are equivalent under rotations and flips of their vertices (Fig- 
ure 3.1), so it suffices to show that fir is invariant under contraction of an 
untwisted edge. This operation is depicted in Figure 3.6 
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Figure 3.5. Computation of /jp for a non-orientable 
Mobius graph T. Here x(Sr) = 0. The surface Sp is a Klein 
bottle. The contributions ±1 add to yield /xr = 4 = (2 — /5) 2 
(for /3 = 4). 
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Figure 3.6. Edge contraction 



Figure 3.7. The diagonal flip operation for a triangulation. 

This relation is obvious since the products of units at either vertex must 
be real and the Wick contraction corresponding to an untwisted edge with 
units e a at either end comes with an overall factor +1. We remark that an 
analogous relation involving a vertex flip for the contraction of a twisted 
edge follows. 

We now compute the invariant /x r for a graph of arbitrary topology. Since 
fir is topological, we may employ a standard graph for each distinct topol- 
ogy labeled by the orientability t| r , number of faces f r and genus g(Sr). 
The key theorem from topology we need here is the connectivity of the 
space of all triangulations of a compact surface with a fixed number of ver- 
tices under the action of a diagonal flip [8] (see Figure 3.7). A diagonal 
flip of a triangulation is exactly the fusion move of the dual Mobius graph. 
For orientable surfaces, the connectivity in concern implies the path con- 
nectivity of the moduli space 9Jt 9>n of smooth complete algebraic curves 
defined over C with n marked points. The connectivity of triangulations of 
non-orientable surfaces has been also established. 

Now note that because of the invariance of /x r under edge contraction, we 
can expand all vertices of valence greater than 3 and contract all vertices of 
valence 1 and 2 to create a trivalent graph. Since the dual of a trivalent 
graph is a triangulation of the surface and the number of vertices of the 
triangulation is the number n = fr of faces of the Mobius graph T, the 
connectivity of the space of triangulations implies that our invariant /x r is 
a constant for each topology of a given n-punctured surface. In the actual 
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computation, it is convenient to use the following representatives for each 
of the three topological classes: 

1. Orientable (tj = 1); a standard graph is given in Figure 3.8. 




Figure 3.8. A standard Mobius graph that represents an 
orientable surface of genus g with n marked points. It con- 
sists of n — 1 tadpoles on the left and g bi-petal flowers on 
the right. 

2. Non-orientable (tj = —1), odd genus; with standard graph given in 
Figure 3.9. 




Figure 3.9. A standard Mobius graph drawn on a non- 
orientable surface of genus g = 2k + 1 with n marked points. 
In addition to n — 1 tadpoles on the left it has k orientable 
bi-petal flowers on the right and a non-orientable one added 
at the top. 

3. Non-orientable (tj = —1), even genus; a standard graph is in Fig- 
ure 3.10. 

Finally, it is easy to calculate /xr for each one particle irreducible compo- 
nents (subgraphs which remain connected after cutting a single edge): 
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Figure 3.10. A standard Mobius graph representing a 
non-orientable surface of genus g — 2k with n marked 
points. It has n — 1 orientable tadpoles on the left, k bi-petal 
orientable flowers on the right and a non-orientable tadpole 
added at the top. 




Each of the above calculations is rather similar: A single line emitted from a 
tadpole subgraph can only ever carry the unit e = 1. As an example, (3.25) 
arises from (i) placing all l's on remaining four lines (ii) l's on one loop 
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and imaginary units on the other which can be done in 2(/3 — 1) ways, 
(iii) the same imaginary unit on both loops giving (3 — 1 possibilities, (iv) 
different imaginary units which can only be achieved for the quaternions 
and in 6 = {(3 — 1){(3 — 2) different ways incurring a minus sign since 
_1 = ijij (say). 

Therefore we find 

(3.28) /ir = (-4 + 6(3 - ^M*-**^) ( 2 - (3)^(3^ . 

This result combined with the first Lemma above completes the proof of 
our main theorem. 



4. Examples 

A simple consistency check is the case f3 = 2. Clearly, non-orientable 
graphs give a vanishing contribution as required for hermitian matrix in- 
tegrals. To obtain a more conventional normalization, we rescale X — > 
f3 l l 2 X in (2. 1) and absorb all but one power of (3 in the couplings t, so that 



(4.1) log 



V 



' fidX] w exp (- | tr X 2 + E°° =1 f tr 
|[^] (/3) exp(-ftr^) 

(_ 4 + 6/3 - ^-^r-hiSr) {2 _ pfrfVh)-!^ 



r r© l Aut ( r )l 



Hence when (3 — 2 
(4.2) 

. / /[^] (2) exp(-ltr^ + E-^trX^) \ 2N fr 

0g { |[^] (2) exp(-|tr^) ) ^|Aut(r)|ll^ 



0) 

r 



We note that 



2A Jr T-r „C?> v iV /r -n- «C?) 



(43) |Aut(r)| 11^/ - 1] | Aut ^ r )| IK r 

re© 1 v 71 i re?H 1 v 71 j 

where y\ denotes the set of all connected ribbon graphs and Aut^r) the 
automorphism group of a ribbon graph disallowing orientation-reversing 
automorphisms. To see (4.3), let T be an oriented ribbon graph. Then either 
(a) T and its flip T are isomorphic as a ribbon graph; or (b) they are different 
ribbon graphs. In case (a), we have |Aut(r)| = 2|Aut< H (T)|. If (b) is the 
case, then Aut(r) = Aut<n(r), but T and f appear as different graphs on 
the right hand side while they are the same as a Mobius graph in <&. 
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Together, equations (4.3) and (4.2) constitute the well-known result for 
the hermitian matrix integral [3]. 

A less trivial test of our new result for the (3 — 4, GSE model is a com- 
parison with the Penner model. To begin with we re-express (3.8) in yet 
another normalization, together with a specialization ti — t 2 — 0: 



(4.4) log 



'j[dX] {4) exp (- \ tr X 2 + Ef=3 I tr xj) 
J[rfX] (4) exp(-ItrX 2 ) 



re© 1 v 71 i>3 



Here we have employed the useful identity 

(4.5) (-l) Sr = (-l) x(Sr) . 
The Penner substitution 

(4.6) tj^-zV 2 - 1 , j>3 
yields the graphical expansion 



(4.7) lim log 



7[^] (4) ex P (-E?= 2 ^ i tr^) 



J[rfX] (4) ex P (-ItrX 2 ) 
(~l) er 

|Aut(r)| 



re© 



Equation (4.7) is valid as an asymptotic expansion of the integral for z — > 
while keeping z > 0, as an element of a formal power series ring (Q[iV])[[z]]. 

The integral on the left hand side can be evaluated explicitly (see (7.2) 
and (7.3) of Section 7). The leading terms as an expansion in z are 



(4.8) lim log 



7[rfX] (4) ex P (-E?: 2 ^trX^) 



J[rfX] (4) ex P (-itrX 2 

= {-T2 N -l N2 + l N3 ) z + °^- 



DUALITY FOR MATRIX INTEGRALS 

On the other hand, our graphical expansion yields 



(4.9) £ 



\A^)\ 



re© 
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+ 




[-2N)z 



(-T2 N -l N2 + l N3 ) z + °^- 



Needless to say agreement is perfect. In fact, as we shall show in Section 7, 
agreement to all orders amounts to known results for the orbifold Euler 
characteristic of the moduli space of real algebraic curves. 

5. DUALITY FOR MATRIX INTEGRALS 

An additional change of variables X — > N l l 2 X in (4.1), absorption of 
all but a single N in the couplings t as well as the substitution 

(5.1) (3 = 2a, 



yields 



(5.2) log 



'l[dX\ 2a) exp (- ^ tr X* + £f =1 ^ tr X*) 



J[dX] 



(2a) 



exp 



x(Sr) 



(3 — a 1 — a) 



re© 



■ ^ tr X 2 

rSr-|x(Sr), 

|Aut(T)| 



-1/2 _ a l/2 



n 



„(j) 
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This formula is invariant under 

(5.3) a — ► a' 1 and N — ► -«JV . 

Remark. 

1. The duality holds graph by graph 3 . 

2. The a — 1 GUE model is self-dual since x(Sr) is even for orientable 
graphs. 

3. The graphical expansion of the a = 2, N x N GSE model is identical 
to that of the a = 1/2 GOE model if the size of the matrices are 
doubled and the contribution of every Mobius graph embedded in a 
non-orientable surface of odd Euler characteristic is multiplied by — 1. 

One might wonder whether matrix integrals exist whose graphical ex- 
pansion coincides exactly with the image of (5.2) under the duality (5.3). 
Although we have no definite answer to this question at this point, in the 
next Section we show that the combination of Poincare duality and the one 
discovered here underly the dualities for correlations of characteristic poly- 
nomials [2, 4, 5, 13]. 



6. Characteristic Polynomial Duality 

The average of products of characteristic polynomials obey dualities be- 
tween GOE and GSE models [5]: 

f[dS\N* N exp f - f tr S 2 ) nj=i det NxN (*i - S) 
(6.1) { -' 



J[d^*"exp(-? trS 2 ) 



/[dX]J* fc exp(-Ntr X 2 ) Hdetf xfe (A -^1) 
J[dX}^ k exp(-N tr X 2 ) 



as well as a self duality for the GUE case 4 [4]: 

J[dX]^ N exp (- f tr X 2 ) nj=i det NxN (*i - X) 



(6.2) 



/ MX]™ exp (- f trX 2 ) 
J[dY]^ exp [- fir Y 2 ) detf xfc (A - y^T Y) 
J[dY]** k expf-f tr Y 2 ) 



3 Physicists would call this a T-duality - valid order by order in perturbation theory. 

4 One might question the reality of the integrals on the right hand sides of (6.2) and (6.1) 
since an explicit \J — 1 appears in the determinants. It is clear, however, from both the 
graphical expansions below and the original derivation in [4, 5] that both integrals are real. 
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The duality relates expectations of /c-fold products of distinct characteristic 
polynomials of N x N matrices to averages over the iVth power of determi- 
nants of certain k x k matrices. Here, A = diag(Ai, ... , A&) is a diagonal 
k x k matrix of real entries. The quaternionic determinant H det in (6.1) is 
defined by 

(6.3) Hdet fcxfc M = detj£ 2 x2fc C(M) , 

where the 2k x 2k matrix C(M) is obtained from the k x k quaternion 
valued matrix M by replacing the quaternionic units by their representation 
in terms of Pauli matrices 1 — > I 2X 2, e.% — » i&i (i = 1, 2, 3). 

To begin with, we demonstrate that the N-k duality for GUE models 
follows from the usual ribbon graph expansion along with Poincare duality 
of graphs on a compact oriented surface. The first step is to represent the 
determinants as vertices of the graphical expansion. Let us assume that the 
parameter X e satisfies \ e > A > for every I and some positive A, and let 
Q\ denote the set of all N x iV hermitian matrices whose eigenvalues are 
contained in the bounded interval [— A, A] . Then for every X e we have 
a convergent power series expansion in A^ 1 : 




where Q c x is the complement of Q\ in the space of all iV x N hermitian 
matrices. Since il\ is a compact space, the first integral on the right hand 
side of (6.4) is a convergent power series in Aj" 1 , A 2 1 , . . . , A fc Set tj = 
— tr A"- 7 . Then Re(^) < 0, and as A — > +oo, tj — > 0. Thus the ribbon 
graph expansion provides each coefficient of the power series expansion of 
this integral in tj as A — > +oo. The second integral on the right hand side 
of (6.4) is a polynomial in A^ 1 whose coefficients converge to as A goes 
to infinity since Q x approaches the empty set and the integrand is bounded. 
Therefore, we obtain an asymptotic expansion formula 
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(6.5) log 



x_ 

A, 



(f[dX]^ N exp (- f tr X 2 ) U.U det NxN (/ - 
\ J[rfX] ( ^exp(-f trX 2 ) 

G(Q[iV])[[Ar 1 ,...,A fc 1 ]] 

The computation of 

7 W ( fc 2 * fe exp (- f tr F 2 ) det^J - ^FA" 1 )' 



(6.6) log 



/[dY]*** exp(-f tr F 2 ) 



can be performed similarly: First we decompose the space of all k x k 
hermitian matrices into two pieces, one consisting of matrices with eigen- 
values in [—A, A], and the other its complement. If A^ > A for every i, then 
det^J — >/— T FA -1 ) can be expanded as before. Asymptotically as an 
element of (Q[X])[[A^ 1 , . . . , \ k % we have 



(6.7) J [dY]^ k exp (- y tr F 2 ) det(J - V^T YA^f 
= J [dY}^ k exp (- j tr F 2 ) exp (^~ N Y1 ^y^ tr (Y^ j 



The appearance of the term tr (FA 1 ) m instead of tr Y m occurring in (6.7) 
replaces products of traces over identity matrices 



/r 0) 



(6.8) N fr = Y[(trl j ) 

j 

incurred in (3.8) as one travels around each face of the graph T, by 

(6.9) JJ(tr A- j ) f r } . 
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(Recall that denotes the number of j-gons in the cell-decomposition of 
S r defined by the graph T.) Therefore, 

(6.10) 

'f[dY]** k exp (- f tr F 2 ) det N (I - y/-[ FA -1 )' 



J[dY]^ k exp (- f tr F 2 ) 
'f[dY]** k exp (- f tr F 2 ) exp (-N £\ ^y^tr (FA _1 )- ? j 
f[dY]** k exp f- f tr F 2 



log 



= log 



where we used 

( — l^r-er _ ^_;Qx(Sr)-/r _ ( — 1)-^ 

Let us denote by T* the dual graph of a ribbon graph T drawn on a compact 
oriented surface S r . We note that Aut^r) = Aut^r*) and 

(6.11) I e r = e r *, 

Since one and two valent vertices are included in the set of ribbon graphs 
d\, the map 

* : m — >9\ 

is a bijection. [Contrast this situation to the Penner model in Section 4, 
where the couplings ti = t 2 = and Poincare duality does not apply.] 
Therefore, 



r€M |Aut OT (r)| 



<°2 



(6.i2) =E^(^(-i)-^---ri(^-r 

This implies that the matrix integrals (6.5) and (6.6) have the same asymp- 
totic expansion. 
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The N-k duality in equation (6.2) is a polynomial identity of degree Nk 
in (Q[iV])[Ai, A2, • • • , Afc], where we define deg(A^) = 1. We must now 
consider also disconnected graphs, since there is no logarithm. The co- 
efficient of the degree Nk — d term of (6.2) is therefore determined by a 
partition d = 2(ei + e 2 H h e m ) corresponding to the product of m con- 
nected graphs consisting of edges. The contributions of connected graphs 
are computed in (6.5) and (6.6). We note that the duality (6.12) holds for 
every surface even when the number of edges is fixed. Therefore, the as- 
ymptotic equality we have derived implies the polynomial identity (6.2). In 
other words, the N-k duality of [4] is a simple consequence of the Poincare 
duality of graphs on a compact oriented surface. 

Our derivation of the characteristic polynomial duality between the GOE 
and GSE models goes quite similarly. Here again we see that the duality is a 
consequence of our graphical expansion formula (3.8) and Poincare duality: 

Using the same trick for characteristic polynomials as in the GUE case, 
from our expansion formula (3.8) we obtain an asymptotic expansion for- 
mula for the GOE side of the duality 

(6.13) 

, 'f[d%? N exp (- f tr S*) Ylti det (/ - 



log 



re© 



f[^]^exp(-f trS* 
'f[dS}^ N exp (- f tr S 2 ) ex P (- £°° =1 ^tr &) 
J[^] ( ^exp(-f trS 2 ) 

, _ , r r9«r-er /v/r-er TTf^r \-3Y' h) 

|Aut(r)| 



(Note the non-standard normalization of the Gaussian exponent yields the 
factor 2~ er .) Its GSE counterpart requires some care: The characteristic 
polynomial of a k x k quaternionic matrix X is defined by 

H det(A -X) = det 1/2 (A/ 2fcx2fe - C{X)) 

and 

trX J = -trC(X) j . 
Thus if all eigenvalues of X are in [—A, A] and A^ > A > 0, then 

( 1 °° 2 

Hdet(7 - XA" 1 ) = exp \ -- ^ - tr (XA"^ 
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Therefore, we have an asymptotic expansion 
(6.14) 

7Wf exp(-iVtrX 2 ) Hdetf xfc (J- v^TXA- 1 )^ 

log 1 



J[dX}^ k exp(-N trX 2 ) 



= log 



'f[dX]^e W (-f tr X 2 ) exp (-N £°° =1 (XA~ 



E 1 i 1 \Sr+er+fro.fr-er JVfVr-er TJ(+ r X-j)/^ 
|Aut(r)| v ; " iil ' 



re© 



£ T4^(- 1 ' /r2/r " riV " r " er n> A "') /?> . 



re .|Aut(r)l 

where we have used the fact that 

that follows from (4.5). (In addition the non-standard Gaussian exponent 
normalization now accounts for the absence of explicit factors k in the 
graphical expansion.) We now see that (6.13) and (6.14) are equal again 
through the dual construction of a Mobius graph. 

The polynomial identity (6.1) follows from the equality of the asymptotic 
expansions. This time each term of (6.1) may have contributions from both 
orientable and non-orientable graphs, but since the dual graph construction 
works for each surface, the equality holds. 



7. The Penner Model 

The Penner model for the hermitian matrix integral provides an effective 
tool to compute the orbifold Euler characteristic of the moduli space of 
smooth algebraic curves defined over C with an arbitrary number of marked 
points [16], [7]. It was discovered in [6] that the Penner model of the real 
symmetric (or GOE) matrix integral yields the orbifold Euler characteristic 
of the moduli spaces of real algebraic curves. Therefore, it is natural to ask 
what topological information is contained in the corresponding Penner type 
model for the quaternionic self-adjoint matrix integral? 

In this Section, we show via our duality (5.3), that the symplectic Pen- 
ner model is identical to the GOE case, except for the matrix size and an 
overall sign for contributions of non-orientable surfaces. The answer to 
the above question is in the negative, the symplectic Penner model offers 
no new topological insight. As we shall see, of the three main classes 
of Mobius graphs -oriented, non-orientable odd x(S , r ) and non-orientable 
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even x(Sr) - only the first two survive the Penner substitution for the cou- 
plings, or in other words, the orbifold Euler characteristic vanishes when 
x(Sr) is even. Therefore the third symplectic Penner type model is not an 
independent topological quantity. 

We also show that the generalized Penner model expressed in terms of 
Vandermonde determinants to powers in 2(N U 1/N) exhibits an extended 
duality that agrees with (5.3) when the power of the Vandermonde is re- 
stricted to 1, 2 or 4. Many explicit formulae and derivations are reserved for 
Appendix B. 

The symplectic Penner model introduced in Section 4 reads 



to be viewed as an element of the formal power series ring (Q[A^])[[^]]. 
This integral is indeed explicitly computable. Symplectic invariance of the 
measure and integrand allows us to diagonalize the matrix variable X — > 

diag(/ci, k 2 , ■ ■ ■ , kx) so that: 



is the Vandermonde determinant. Using the asymptotic expansion tech- 
nique established in [14], the Selberg integral formula and the Stirling for- 
mula an explicit asymptotic expansion, even valid for every a E N, can be 
computed for the integral 




(7.2) 




where 



i<j 



(7.3) 



K(z, N, a) = lim log 



( 



^A 2 "WIir = iexp(-E. 





J RJ vA 2a Wnf =1 exp(-f) 
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(In what follows K(z, N, a) should be regarded as as the asymptotic ex- 
pansion, not the underlying integral; an explicit expression is given in Ap- 
pendix B.) Specializing to the a = 1 hermitian case yields a very compact 
result corresponding to the original formula of Penner [16]: 

K{2N1)= E e^aiV f.,)^ 

(7.4) 2-2g-n<0 



E 



(—lY r 



, fe9t |Aut*(r)| 



Identifying ra with / r and g as the genus, yields the well known generat- 
ing function of the Euler characteristic x(%Rg,n) of the moduli of complex 
algebraic curves of genus g and n marked points. 

For the a = 2 symplectic case a similar simplification occurs and (7.3) 
can be written in terms of the a — 1 hermitian result plus additional terms 
corresponding to non-orientable surfaces of even genus g — 2q with m + 
1 — 2q marked points 



(7.5) 

K(z,N,2) 



K(z,2N,l) 

(2q + n 



E 



2)\{2 2q ~ l - 1) 



(?>0,n>0 
l-2<?-n<0 



(2q)\n\ 



b 2q (2N) n (-z) 



n(_ \2q+n-l 



Comparing (7.5) with (7.1), we obtain 



(-l) er 1 (2g + n-2)!(2 2 ^-l) 
1 } ^ |Aut(r)| 2 (2q)\n\ 2? ' 

r: / r =n, <y(5 r )=2g 1 V 71 V yy 

where the summation is over all connected non-orientable Mobius graphs 
with n > faces that are drawn on a non-orientable surface of genus 2q 
satisfying a hyperbolicity condition l — 2q — n< 0. The formula (7.6) is in 
exact agreement with the formula for the orbifold Euler characteristic of the 
moduli space of smooth real algebraic curves of genus 2q with n marked 
points that can be found in [6] and [15]. 

To study the GOE-GSE duality for the Penner model, we need the ex- 
pansion of the analog of formula (7.3) valid for the single power of the 
Vandermonde determinant relevant to the GOE model. Indeed an integral 



28 



M. MULASE AND A. WALDRON 



formula for 



(7.7) 



J(z, N, 7) = lim log 



V J R ,v|A(A:)p/7nf =1 exp(-f)^ 



valid for every positive integer 7 e N (z'.e. for all powers 2/N of the Van- 
dermonde) was derived in [6] in order to compute (7.6). (See Appendix B.) 

The matrix integrals K{z, N, a) and J(z, N, 7) are closely related: Ob- 
viously 



expresses the 7 = 2 GOE case in terms of the 7 = 1 oriented hermitian 
result and a sum over non-orientable contributions with odd Euler char- 
acteristic. Notice, as claimed above, in the graphical expansions of the 
orthogonal Penner model J(2z, 2N, 2) and the symplectic Penner model 
K{z, N, 2), non-orientable surfaces of odd genera do not contribute. This 
corresponds to the fact that the orbifold Euler characteristic of the moduli 
space of smooth real algebraic curves of odd genus 2q + 1 with n marked 
points is for any value of q > and n > 0. 

More importantly, observe that the GSE and GOE formulae (7.5) and (7.9) 
almost coincide except that the GOE formula is for matrix size 2N and the 
non-orientable odd Euler characteristic terms differ by an overall sign. (The 
appearance of 2z rather than z will be cured by the appropriate normaliza- 
tion given below and in the master formula (5.2).) This is precisely our 
duality (5.3). 

Finally, we show that the duality between GOE and GSE extends to ar- 
bitrary positive integers for the two types of Penner integrals introduced in 
this Section. Let r e NUl/N and set 



(7.8) 



J(z,N,l) = K(z,N,l). 



The identity 



(7.9) 




(7.10) I(z,N,r) 
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= Hill log -; -yr 

m ^°° V U\ A ( k )\ 2r Ul 1 ^p{- k t)dk l 

G(Q[iV,iV- 1 ,r,r- 1 ])[[z]]. 

Then we have 

( K(^,N,a) , r = a G N , 

(7.11) J(^,iV,r) = ^ 

[ J(f,iV, 7 ) , r = r 1 Gl/N. 

From inspection of the explicit asymptotic expansion formulae of (B.l) 
and (B.9) presented in Appendix B, we obtain an extended duality 

(7.12) I(z,N,r) = 7(z, -rTV, r -1 ) 

for an arbitrary positive integer r. This is in agreement with our duality (5.3) 
for r = 1,2. 

8. The Central Limit Theorem 

To prove a central limit theorem for large matrix size N, we need to 
show that the leading dependence is Gaussian in the coupling constants tj. 
More precisely, define the Gaussian expectation value of f(X) for GOE 
(a = 1/2), GUE (a = 1), and GSE (a = 2) as 

J[dX] {2a) exp(-^trX 2 ) /(X) 
(8.1) / (X Wa) = — ™ ^— ^— , 

and consider 

; e xp(E,^trX^) 



(8.2) V(t,N,a) = 



(N,a) 



exp (E, ^ (tr Xi 



(N,a) J 

The expansion formula (5.2) shows that the contribution of a connected 
Mobius graph r G <5 to log /exp f V, ^ tr X*') \ is 

\ V J J / / (N,a) 

(8.3) 

|Aut(r)| l/' • 

Since 
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is the sum of all contributions from 1 -vertex Mobius graphs, we see that 
logV(t, N, a) has no terms coming from 1-vertex graphs. In particular, it 
has no terms with a positive power of N. Indeed, the power of N in (8.3) 
is strictly positive only when x(Sr) = 2 and v r = 1, i.e. T is an orientable 
planar graph with one vertex. Therefore, limjv^oo log V(t, N, a) consists of 
contributions from graphs that have two or more vertices and x(Sr) — v r = 
0. But this is possible only when x(Sr) — fr — 2. In other words, T is an 
orientable planar graph with exactly two vertices contributing (see (8.3)) 

2a 

|Aut(r)|^ 2 ' 

where j 1 and j 2 are the valences of the two vertices of T. Altogether, we 
have established the following 

Theorem 8.1. The Central Limit Theorem for GOE (a = 1/2), GUE (a = 
1) and GSE (a — 2) ensembles is 

[OX]) | / , — Tl' a J | , 

(N,a) 



lim log 



(8.4) 



E 



a 



lAuwni 

T connected, oriented, planar 
2-vertex ribbon graph 

where j\ and we the valences of the two vertices of the ribbon graph Y. 

We notice that the formula is the same for all three ensembles except for 
the overall factor of a. In particular, only oriented planar ribbon graphs 
contribute in the large limit. This mechanism was observed long ago by 
't Hooft in the hope that large A^ quantum chromodynamics could be solved 
exactly [19]. 

9. Conclusions 

Let us tabulate the patina of results gathered here: 

• The asymptotic expansion of the three Gaussian random matrix en- 
sembles is expressed as a sum over Mobius graphs. 

• These expansions are related by a duality 

(9.1) a — > a; -1 and AT — ► -aN 

The a — 1 GUE model is self-dual and sums over only ribbon graphs. 
The duality between a — 1/2 GOE and a = 2 GSE models amounts 
to an equality of graphical expansions up to a factor 

(_l)x(Sr) 
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for any graph T. 

• When specialized to Penner model couplings, the Selberg integral rep- 
resentation yields an asymptotic expansion for all a e N U 1/N and 
the duality (9.1) holds for this extended set of a's. 

Therefore, the first and probably most interesting question one might pose, 
is whether our graphical expansion formulae can also generalized to the ex- 
tended set of a E N U 1/N. I.e.; 




Let us briefly postpone a discussion of this issue while enumerating several 
other questions for which we have no immediate answers: 

1. Do there exist matrix models where the duality holds exactly, without 
a factor (-1)*^? 

2. What is the significance of the minus sign in the transformation 

N -> -aN? 

Is there an interpretation where traces over N x N matrices are re- 
placed by a supertrace and in turn bosonic matrix integrals by fermionic 
ones? 

3. Why the factor a in the transformation 

N -> -aN? 

Is there a generalization of the GSE dual to the GOE at odd values of 

After this disquisitive interlude, we return to the postponed question. Let 
us examine the generality of Lemma 3.3 which claimed that a topologi- 
cal invariant of a punctured surface S r was obtained by counting signed 
configurations of units {l,ei . . . , e 2a -i} on the associated graph T. Its 
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proof relied on the following: (i) The units all square to ±1 and {±1, ie^} 
is a group, (ii) At any vertex, their product was ±1 and therefore real, 
(iii) Units whose square was (— )1 were joined by Wick contractions of 
(anti)symmetric matrices. 

Therefore generalized matrix models whose graphical expansion is com- 
pletely determined via our methods can be written down based on a larger 
group of "imaginary" units {±l,±/ a ;±ej : /„ = l,ej = — 1}. Sim- 
ple examples are generated by considering elements drawn from Clifford 
algebras. I.e., the Pauli matrix representation of the quaternions can be gen- 
eralized to larger sets of higher dimensional Dirac matrices. Although these 
Clifford type models seem not to generate theories with the new values of a 
exhibited in the Penner model, it would be interesting to investigate whether 
new matrix models of this type can indeed be constructed. 
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Appendix A. Quaternionic Feynman Calculus 

Since the quaternions are the last real associative division algebra, it is 
natural to develop a manifestly quaternionic Feynman calculus. Again, let 
us consider GOE, GUE and GSE models all at once via the unified notation 

p-i 

(A.l) X = S + Y,e i A i 

i=i 

where the x A r matrices X are real, complex or quaternionic self-adjoint 

(A.2) Xt EE X T = X , 6j EE 6j , 

depending on the value of f3 — 1, 2 or 4, respectively. 
To begin, we need a shift identity 

(A.3) expQtr(g T X + Xg T )^) f(Q) = f(X + Q), 

where the "quantum variable" 

(A.4) Q = S + ^e i A i = Q* 
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and the N x N matrix of derivatives d is given by 

i / 

a = b . 



(A.5) 



d, 



oh 



2 \dS ab 
d 



Q S aa 



The identity (A. 3) holds thanks to the commutation relations 



(A.6) 

Note also that 
(A.7) 



itr (<9 T X + X«9 T ),Q 



X. 



d,^ti(Q T X + XQ T )] = Q, 



although [d, Q] = I only for the real and complex cases (3 = 1,2. 
We may now rewrite matrix integration as differentiation 5 

[dX] (4) exp (- - tr X 2 + J2 j tr Xj ) 

oo 

[rfX] (4) exp (- - tr X 2 + £ | tr (X + Qf) 



3=1 



Q=0 



= J [dX] (4) exp(- itr (x - «9 T ) 2 + 1 tr «9 t2 ) exp(^ *j tr Q 
= (/[^] (4 fxp(- 



['•' X i , , x n [ - - tr X 2 ) exp (- tr d 2 ) exp ( £ 4 tr Q J 

' .7 = 1 ^ 



Q=0 



Q=0 



The first factor on the last line is just an overall normalization while the 
two exponentials can be expanded in terms of Feynman diagrams: the n-th 
order term in the expansion of each exponential is interpreted as either n 
edges or n vertices, respectively. 

Let us give some details: The operator | tr d 2 acting on a quantum vari- 
able, yields 



(A.9) 



[- trd 2 ,Q}=d 1 



which is represented graphically as attaching a ribbon edge to a vertex since 
the operator | tr d 2 is to be viewed as an edge. Note that by this rule, a 



5 The result is equivalent to the one obtained using a Schwinger source term. The refor- 
mulation presented here is often called the background field formalism, a simple account 
may be found in the on-line textbook [17]. 
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vertex emitting a Q ab is replaced with one emitting d ba which amounts to a 
twist. Attaching the other end to an adjacent vertex yields 

(A. 10) (d ab Q cd ^ = ±p 5 ac 5 bd + 1(2-13) 5 ad 5 bc . 

Note that the brackets on the left hand sides above indicate that we are 
computing the derivative rather allowing it to continue acting to the right 
as an operator. In particular, observe that for the GUE case, (3 = 2 so no 
twisted ribbon graphs can appear. 

For the GOE and GUE models we are done, one simply attaches all possi- 
ble edges to vertices using the above rules and finds the usual known results. 
The symplectic case is more subtle however, thanks to the quaternionic non- 
commutativity of d and Q. In particular 

(A.ll) (d ab f(Q) Q cd ) + (d ab f{Q)) Q cd + f(Q) (d ab Q cd ) . 

for some function / of the quaternionic matrix Q. I.e., the quaternion valued 
operator d does not satisfy the Leibniz rule. However, a generalized Leibniz 
rule does apply: First note that for any 

(A. 12) Q = S + Y,e a A a 

a 

where the real matrices S and A a need not have any definite symmetry 
properties (so that Q is not necessarily self adjoint) we have 6 

(A.13) ^2 e aQ e a = — Q + (2 — (3) ~Q . 

a 

Therefore we have the generalized Leibniz rule 
(A. 14) 

dab(Q Qcd} - (dabQj Qcd = S ac S bd Q - -(2 - (3)(5 ac 5 bd - 5 ad 5 bc )Q . 

Note that for (3 = 1,2 the right hand side is equal to Q (d ab Q c d} , expressing 
the commutativity of real and complex numbers. This relation is the central 
graphical rule for our quaternionic Feynman calculus and is depicted in 
Figure A.l. Note that it reverts to the rule (A. 10) when Q = 1. 

It is now possible to compute any term in the expansion of (A. 8) in terms 
of graphs. However, for quaternionic matrices, when connecting vertices 
with ribbon edges, intermediate vertices and unconnected edges may be 
twisted and/or flipped according to (A.l). It is easy, but tedious to verify 
that the results for simple graphs coincide with our general formula (5.2). 



6 When (3 = 1, the sum on the left hand side is empty and equal to zero, while the right 
hand side vanishes for real Q. 
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Figure A.l. Quaternionic Feynman rule. Observe how 
connecting vertices with edges changes vertices yet to be 
connected. 

Appendix B. Generalized Penner Model 
In this Appendix, we derive the asymptotic expansion formula 

(B.l) 



Ln A 2a (fc) []£=! exp - 



K(z, N, a) = lim log 

m— *oo 

h ± 1 

— nz 2 ™- 1 + y - , 

- 1) ^ 4m 
-a m (a'-' 2q - 1)N 



E 



^ 2m (2m 



/ R ,A^(fc) nil exp 

l) m a m N m z m 



2 A (2g)!(m+ 1 - 2g)l 



m=l ij=0 

r m ] r m+1 i 



oo L 2 - 



EE E 

m=l 9=0 s=0 



-l) m (m - l)\b 2q b 



2.s 



(2g)!(2s)!(m + 2-2g-2s)! 



2<jjym+2— 2g— 2s^?ti 



which is valid for every positive integer a. Here the 6 n 's are the Bernoulli 
numbers defined by 
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The key techniques are the Selberg integration formula, Stirling's formula 
for T(l/z) and the asymptotic analysis of [14]. First we note that as an 
asymptotic series in z when z — > while keeping z > 0, we have 



lim log | / A 2a (k) ff exp ( - V k{]dki 

= log | f^eiz^z 2 ^ / A 2Q (A;)TTe fei A;/dA; i 



'[0,oo)' 

(For the mechanism changing the integration from [R^ to [0, oo)^, we refer 
to [14].) This integral can be calculated by the Selberg integration formula: 



/ 



N N-l 



[o,oo)" f=l t=o T(l + a) 



Therefore, we have 

/ /• AT 2m j/2-1 

(B.3) lim log / A 2a (k) TT exp ( - V ) rfj^ 

— yj RN a. i v ^ j y 

AT, N N, aJV(JV-l), 
= c + — log z H 1 log z H log z 

2 z z 2 

AT-1 



+ \ogl[r(l + -+jaY 

3=0 ^ " ' 



where c is the constant term independent of z. Since (B.l) does not have 
any constant term relative to z, here and below we ignore all constant terms 
independent of z (but possibly iV dependent). The product of T-functions 
can be calculated by the recursion formula, noticing that a is an integer: 

N-l , 
3=0 V 

AT-1 ia 



(B.4) 



i=0 j=0 V J 
/, \ AT AT-1 a , . Af-i-1 

n^riS) 1 1 II ( ' /) 

w 0)" nil ( 1+2(1+( ;- 1)Q+J) ) 



i=0 j=l 

jVAT-la-1 /, , /-,,/• , \ .\\ N-i 



i=l j=0 
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We now apply Stirling's formula for logr(l/z) to obtain, up to a constant 
term: 

N-l 



log Yl rfl + -+ja) 

3=0 ^ ' 



N N. ^ b 2m 



= \ogz + — log z + V - 2m 

rx z z 2 ^2m2m-l 

(B.5) m=1 

- iV log z log z 



\m m 



+ E E E ^(- 1 ) m ~ 1 ( iV - oa + (< - 1)« + i) 
lit 

m=l i=l j=0 

We note that all negative powers of z and log z related terms in (B.3) cancel 
out using (B.5). Finally, we obtain 

(B.6) K ( z,N, a) = Y.^t^Tf^ 

m=l 

oo N—l a 1 

+ E E E -(-ir^iN -i-i)( ia + 3 ) m z m . 

TYl 

m=l i=0 j=l 

This last sum of powers can be calculated using Bernoulli polynomials, 
from which (B.l) follows. 

Using a formula for Bernoulli numbers 



n /2n \ n /2n\ 

(1 - 2n)b 2n = {£)b 2q b 2n - 2q = [ 2 JS^- 2 ,2 29 , 



q=0 x " / q=0 

and noting that b 2 = 1/6, we recover the original formula of Penner for 
a = l[16]: 

K(z,N,l) = -Y b -^Nz 2m - 1 
^ 2m 



m=l 



(B.7) 



1 ff (m-l)!(2g-l) m+2 _ 2gr )m 

^ (20 + 71-3)1(20-1) ^ x2 fl +n-2 



2-2g-n<0 

EC — 1 Y r 
1 ; j\r /r C-z) er ~ t ' r 
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For a = 2, (B.l) simplifies again: 
(B.8) 

K(z,N,2) = -Y b -^Nz 2m ~ l 
y ' ' ^ 2m 

m=l 

, 1 \^ \ (m-l)!(2g-l) 9Tra+2 -2g A r m+2-2g/_ ,p 

+ 2 ^(2g)!(m + 2-2g)!^ 2 " ( ^ 

[221 

1 °° 1 2 J / 1\|/q2(J— 1 1\ 

"2^ (2,)!( m+ l-2 9 )! ,,2 « 2 W <- 2) ■ 

Note that the first two lines of (B.8) are identical to the Penner model 

\K{z,2N,l). 

The following integral formula, again valid for every positive integer 7 e 
N, has been established in [6]: 

(B.9) 




| R ,|A(A;)|V7 n f =i exp(-f)^ 



— , / \ m 

ZLm 



... -1 m iV m - 
2m (2m — 1) 7 \7/ 4m \7 

= 1 v ' x 7 m=l 



1 ^ (-1)>-1)S A _ _M ^h-i-2, (A 

2 ^^(2g)!(m+l-2g)! V 7 1 " 2 V W 



^ [22] [ m+1 ] g 

f f V- (-l) m (m-l)\b 2q b 2s _J_ N m+2-2q-2s ( 

^ (2g)!(2 S )!(m + 2-2g-2 S )! 7 1 " 28 V7 

Our duality (7.12) for the generalized Penner model follows from compar- 
ing (B.l) with (B.9). 
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